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Abstract 

ALADIN multifragmentation data show features characteristic of 
a critical behavior, which are very well reproduced by a bond perco- 
lation model. This suggests, in the context of the lattice gas model, 
that fragments are formed at nearly normal nuclear densities and tem- 
peratures corresponding to the Kertesz line. Calculations performed 
with a lattice gas model have shown that similarly good reproduction 
of the data can also be achieved at lower densities, particularly in the 
liquid-gas coexistence region. 

In nuclear multifragmentation studies the presence of a phase transition 
is often deduced, however, its nature is not unambiguously identified. In 
small systems the asymptotic behavior is strongly modified by finite size ef- 
fects so that the distinction between first- and second-order phase transitions 
becomes very difficult. Simulations with one-component lattice gas models 
have shown that critical-like features are observed in finite systems not only 
along the Kertesz line but also inside the liquid-gas coexistence zone, i.e. 
the first order phase transition can mimic critical behavior [I1E1I3]. 

We have performed similar criticality analysis with a more realistic 
isospin dependent lattice gas model (LGM) which includes Coulomb in- 
teractions. Criticality is deduced from fluctuations of the largest fragment 
size. The fluctuations are measured by cumulants (cumulant ratios) of the 
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probability distribution of the largest fragment size, such as the skewness, 
Ks, and the kurtosis excess, K4. As it was shown in percolation studies, the 
pseudocritical point in finite systems is indicated by K3 = and a minimum 
value of K4 [4J. Our LGM results are confronted with the ALADIN S114 
data on fragmentation of ^^'^Au projectiles, following the procedure used in 
percolation analysis [5j. 

In LGM a system is composed of Zq protons and Nq neutrons which 
are distributed on a simple cubic lattice with L sites. The configuration is 
described by the site occupation numbers nj. If the i-th. site is occupied by 
proton (neutron) then rij = 1 (-1), if the site is unoccupied then = 0. 
The system density is given as p = pqAq/L, where pQ is the normal nuclear 
density and Aq = Zq + Nq. The Hamiltonian has the following form: 



i=l <i,j> ni=nj=l, j^i 



where pi is the nucleon momentum, m is the nucleon mass, e^^n is the 
nuclear coupling constant (en = e_i_i = 0, ei_i = e_ii = 5.33 MeV), 
/c = 1.44 MeV fm is the Coulomb constant, and is the distance be- 
tween sites i and j calculated for the lattice spacing tq = 1.8 fm [6l |7j. The 
nuclear interaction acts between nearest neighbors < i,j >. The equilib- 
rium configurations were sampled according to the potential energy with the 
Metropolis algorithm for the canonical ensemble. Given a configuration, the 
momenta of nucleons were randomly drawn from the Maxwell-Boltzmann 
distribution. Clusters were defined using the Pan-Das Gupta prescription: 
two neighboring nucleons belong to the same cluster if their relative kinetic 
energy is insufficient to overcome the attractive bond [8]. Free boundary 
conditions were applied to account for the presence of a surface. Typically, 
a sample of 10^ events was obtained for each considered system at a given 
temperature. 

Let us start the examination of the largest fragment size (charge) distri- 
butions in the vicinity of the thermal critical point. This point is located 
at the density p = O.Spo- The critical temperature is indicated by a max- 
imum of the specific heat Cy. As can be seen from Fig. 1 (the lower left 
panel) the critical temperature is about 4.8 MeV for the investigated sys- 
tems with Zq = 36 -7- 68. In the present calculations we have assumed the 
ratio Aq/Zq = 2.5 of the ^^^Au nucleus. At the critical point, the cumulants 
K3 and K4 exhibit = and minimum of K4 as in the case of the perco- 
lation transition. The presence of such criticality features is also observed 
for systems at other densities. We have performed calculations for differ- 
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Figure 1: Lattice gas model predictions. Left panel: the cuniulants of the Z^ax 
distribution and the specific heat Cy as a function of the temperature at the density 
p = 0.5/90 for different system sizes Zq. Right panel: locations of the critical- like 
points (full circles) in the LGM phase diagram. The liquid-gas coexistence and the 
Kertesz lines are schematically drawn on the basis of [9]. 



ent density values within the range (0.3 ^ 1.0)po- The temperature at the 
critical-like point as a function of the system density is shown in the right 
panel of Fig. 1 for a system with Zq = 72. This result corroborates that 
in finite systems a critical-like behavior may be observed not only along the 
Kertesz line but also at subcritical densities and temperatures. 

In searching for the criticality signal, the control parameter T may be 
replaced by the well measured quantity Zbound- The signal {K^ = and 
minimum K4 of about -1) remains well preserved and appears at a Zbound 
value corresponding to the critical temperature [4J. Figure 2 shows the 
signal locations on the Zfjo^nd axis for percolation, LGM and the ALADIN 
S114 data on fragmentation of ^^''Au projectile spectators at the incident 
energies of 600 -i- 1000 MeV per nucleon |10j . The model calculations were 
performed for different system sizes and the obtained results are plotted 
against Zbound normalized to the system charge Zq. In such a representation 
the signal locations are nearly independent of the system size. In particular. 




Figure 2: The cumulants of the Z^ax distribution as a function of Z^ound- Perco- 
lation and LGM calculations are performed for different assumed system sizes Zq. 
The LGM results are for the density p = Q.5po. The experimental data concern the 
fragmentation of ^^''Au projectiles at relativistic energies. 



the percolation transition is observed at -^bowndZ-^o = 0.84. In the case of 
the experimental data the size of the fragmenting system is not precisely 
established. With decreasing .^bound the average system size systematically 
decreases while the excitation energy per nucleon increases. As indicated 
by the cumulants, the criticality conditions are reached at Zj^ound = 54. To 
observe the criticality signal in percolation at the same position Zhound = 
54, the system charge must be assumed as Zq = 54/0.84 ~ 64. In LGM 
the signal matches the experimental position when Zq ~ 54/0.67 ~ 80 in 
the case p = 0.5po- The LGM system charges for this and other densities 
are shown in Fig. 3. They are evidently too large, since they exceed the 
projectile charge. The percolation result is in a good agreement with an 
experimental estimate. 

In LGM the system sizes are too large with respect to the experimental 
data, due to overpredicted multiplicities of Z = 1 fragments. This does 
not affect Zhound values. Taking into account only fragments with Z > 1, 
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Figure 3: The system size Zq at the critical-hke point located at Zhound = 54. The 
experimental mean value of Zq is estimated as 64.5 ± 4.5 [11]. The full circles show 
the LGM results obtained for different densities. The percolation result is indicated 
by the triangle. 



the fragment charge partitions are in agreement with the data, which can 
be concluded from Fig. 4. The comparisons shown in Fig. 4 are made at 
Zhound = 54 for various observables related to the fragment charge parti- 
tions. The first column recalls the resemblance between the data and the 
percolation results [5|. In the next columns the same data are compared 
with the LGM model predictions for three different densities. In all the 
cases the reproduction of the experimental data is astonishingly good. 

In summary, the largest fragment charge distributions characteristic of 
critical behavior, indicated by = and minimum X4, are observed in 
LGM along the Kertesz line with an extension to subcritical densities and 
temperatures. In the ALADIN data the criticality signal is observed at 
Zhound = 54. The system sizes Zq required by LGM to reproduce the ex- 
perimental signal are significantly larger than the experimental estimate. 
Considering only fragments with Z > 1, the simulated characteristics of 
the fragment charge partitions at the critical-like point are nearly identical 
for different freeze-out densities, and are in a good agreement with the ex- 
perimental data. This indicates that the observed differences in Zq values 
are generated by different multiplicities oi Z = 1 fragments. Despite this 
disagreement, one can conclude that observables related to the fragment 
charge partitions, in particular the largest fragment charge distribution, are 




Figure 4: Comparisons between the ALADIN data (circles) and models predictions 
at Zhound = 54. The lines represent calciilations of the bond percolation inodcil 
(first column) and the lattice gas model for densities indicated in the figure (next 
three columns). Prom top to bottom: the distribution of largest fragment charge, 
the mean fragment multiplicity as a function of the fragment charge (largest frag- 
ment excluded), the mean fragment charge as a function of the fragment rank, the 
probability distribution of the second largest fragment charge, and the multiplicity 
distribution of fragments with Z > 2. 



useful in revealing a phase transition (critical behavior) in multifragmenta- 
tion, however, its location in the phase diagram cannot be unambiguously 
determined. 

Further investigations with LGM are necessary to evaluate how the mul- 
tiplicity of Z = 1 fragments is sensitive to model parameters and cluster 



definition options. It will be also interesting to examine isotopic distribu- 
tions, which might resolve the freeze-out density question. 

This work has been supported by the Polish Ministry of Science and 
Higher Education grant N202 160 32/4308 (2007-2009). 

References 

[1] GULMINELLI F. AND CHOMAZ P., Phys. Rev. Lett, 82 (1999) 1402. 
[2] CARMONA J., RICHERT J. AND WAGNER P., Phys. Lett, B531 
(2002) 71. 

[3] GULMINELLI F. AND CHOMAZ P., Phys. Rev., C71 (2005) 054607. 

[4] BRZYCHCZYK J., Phys. Rev., C73 (2006) 024601. 

[5] BRZYCHCZYK J. ET AL., Proceedings of the International Workshop 

on Multifragmentation and Related Topics (IWM2009), Catania, Raly, 

November 2009 ; larXiv: 1003.23271 [nucl-ex] . 
[6] SAMADDAR S. K. AND GUPTA S. D., Phys. Rev., C61 (2000) 

034610. 

[7] LEHAUT G., GULMINELLI F., AND LOPEZ O., larXiv:0809.2T63] 
[nucl-th] . 

[8] PAN J., GUPTA S. D. AND GRANT M., Phys. Rev. Lett., 80 (1998) 
1182. 

[9] CAMPI X. AND KRIVINE H., Nucl. Phys., A620 (1997) 46. 
[10] SCHUTTAUF A. ET AL., Nucl. Phys., A607 (1996) 457. 
[11] POCHODZALLA J. ET AL., Phys. Rev. Lett, 75 (1995) 1040. 



